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The temperature dependence of the sigma meson and pion masses is studied in the 
framework of the 0{N) model. The Cornwall-Jackiw-Tomboulis formalism is applied 
to derive gap equations for the masses in the Hartree and large- A'^ approximations. 
Renormalization of the gap equations is carried out within the cut-off and counter- 
term renormalization schemes. A consistent renormalization of the gap equations 
within the cut-off scheme is found to be possible only in the large- A'' approximation 
and for a finite value of the cut-off. On the other hand, the counter-term scheme 
allows for a consistent renormalization of both the large- A'' and Hartree approxima- 
tions. In these approximations, the meson masses at a given nonzero temperature 
depend in general on the choice of the cut-off or renormalization scale. As an appli- 
cation, we also discuss the in-medium on-shell decay widths for sigma mesons and 
pious at rest. 

I. INTRODUCTION 

Although chiral symmetry is manifest in the Lagrangian of quantum chromodynamics (QCD) for 
vanishing quark masses, quantum effects break this symmetry spontaneously in the QCD vacuum. 
At temperatures of order 150 MeV, however, lattice QCD results indicate that chiral symmetry is 
restored 0]. Such temperatures are expected to be reached in ultrarelativistic heavy-ion collisions 
at CERN-SPS, BNL-RHIC, and CERN-LHC energies g. The restoration of chiral symmetry may 
lead to observable consequences, for instance, in the dilepton mass spectrum or the formation of 
disoriented chiral condensates 

QCD with Nf massless quark flavors has a SU{Nf)L x SU{Nf)ii symmetry. The order parameter 
for the chiral transition is therefore x'-' = (Oi^jj), i,j = l,..., Nf. For Nf — 2, SU{2)l x 5[/(2)fl is 
isomorphic to 0(4). Consequently, the effective Lagrangian for the order parameter x^'' falls in the 
same universality class as the 0(4) model, with order parameter 0% ? = 1, . . . , 4. Thus, if the chiral 
transition is second order in QCD, the dynamics (and the critical exponents) will be the same as in 
the 0(4) model, provided one is sufficiently close to the transition temperature This motivates 
the study of the 0(4) model |^ as an effective low-energy model for QCD. The study which will be 
presented here is based on the 0{N) model for arbitrary N. To make contact with QCD, however, 
we take = 4 in all numerical calculations. 

At finite temperature, the naive perturbative expansion in powers of the coupling constant breaks 
down, requiring resummation schemes to obtain reliable results Typically, these schemes aim 

to include thermal fluctuations to all orders in the calculation of physical quantities. Over the years, 
many ways to achieve this have been pursued (some approaches are more rigorous, some are more 



or less ad hoc; Refs. [p^ constitute an incomplete list). The present study of the 0{N) model 
focuses on the Hartree approximation and its large- TV limit (the large- TV approximation). 

The Hartree approximation is known from many-body theory and generically represents the self- 
consistent resummation of tadpole diagrams |p^ . There is, however, no unique prescription to 
perform this resummation. This has led to the confusing situation that the same term "Hartree 
approximation" has been used for resummation schemes which actually differ in detail ]lO| , [l^ , p^ . 
In this work, the so-called Cornwall-Jackiw-Tomboulis (CJT) formalism is applied to derive a 
Hartree approximation [20[ , and we shall use the term "Hartree approximation" exclusively for 
the resummation of tadpoles within the CJT formalism. 

The CJT formalism can be viewed as a prescription for computing the effective action of a given 
theory. In general, the CJT formalism resums one-particle irreducible diagrams to all orders. The 
stationarity conditions for the effective action are nothing but Schwinger-Dyson equations for the 
one- and two-point Green's functions of the theory. What is usually flj,^ referred to as the 
Hartree approximation in the context of the CJT formalism is the special case where only one- 
particle irreducible tadpole diagrams are included in the resummation. (To be precise, the original 
work of Cornwall, Jackiw, and Tomboulis referred to this as the "Hartree- Focfc approximation" . ) 
In this case, the equations for the two-point Green's functions simplify to self-consistency conditions, 
or "gap" equations, for the resummed masses of the quasi-particle excitations, i.e., in our case, the 
in-medium sigma and pion masses. 

The 0{N) model has been previously studied using the CJT formalism by Amelino-Camelia [ p^ , 
and Petropoulos . The former work addressed renormalization using the cut-off renormalization 
scheme, but did not present solutions of the gap equations. In the latter work, the gap equations 
were numerically solved, but the issue of renormalization was not treated. In this paper, we complete 
these investigations by discussing the renormalization of the gap equations in the cut-off as well as 
the counter-term renormalization scheme and by presenting the corresponding numerical solutions. 

Renormalization of expressions obtained in self-consistent approximation schemes is non-trivial. In 
perturbation theory, it is sufficient to perform renormalization in the vacuum, T = 0, order by order 
in the coupling constant, as a finite temperature does not introduce new ultraviolet singularities 
p2| . Consequently, the perturbative renormalization of the 0{N) model is straightforward and 
has been known for a long time p^ . Self-consistent approximation schemes, however, in general 
resum only certain classes of diagrams. This has the consequence that performing renormalization 
after resummation may require renormalization constants that are no longer independent of the 
properties of the medium, see Ref. [|lO| and below. A resummation scheme which circumvents this 
problem by renormalizing prior to resummation is the so-called optimized perturbation theory [ pT| . 
This approach will not be discussed here. 

Our main results are the following. In the cut-off renormalization scheme, we find that taking the 
cut-off, A, to infinity the masses of the sigma meson and the pions become identical |p^ , even in 
the phase where chiral symmetry is broken. This is clearly unphysical, as the pions are Goldstone 
bosons and thus much lighter than the sigma meson. Moreover, it is a well-known fact that the 
0{N) model becomes trivial in the limit A ^ oo. Consequently, renormalization of the 0{N) model 
within the cut-off scheme can only be meaningfully studied for a finite value of A. Even then, fixing 
A to give the observed meson masses in the vacuum, we find that in the absence of explicit chiral 
symmetry breaking, the Hartree approximation requires A to vanish. The large- approximation 
does not have this problem and allows for nonzero values of A. 

In the counter-term scheme, renormalization can be performed within both the Hartree as well as 
the large- A'^ approximation. In the Hartree approximation, the value of the renormalization scale, /i, 
is uniquely fixed by the vacuum mass of the sigma meson in the absence of explicit chiral symmetry 
breaking. In the large-A'^ approximation, this constraint does not exist, and can be chosen arbi- 
trarily. As A or /X are free parameters in the large- A'^ approximation, at any given temperature the 
values of the meson masses depend on the choice of these parameters. This is a consequence of the 
fact mentioned above that, for self-consistent resummation schemes, the renormalization constants 
may depend on the temperature. 

The outline of the paper is as follows. In Section II, we briefly discuss the effective potential 
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within the CJT formalisin |2^. In Section III, this formalism is apphed to derive the effective 
potential for the 0{N) model in the Hartree approximation. Section IV is devoted to a discussion 
of the stationarity conditions for the effective potential, which lead to gap equations for the sigma 
and pion masses. The renormalization of the gap equations is then performed in Section V within 
the cut-off and the counter-term schemes. In Section VI we present numerical results. Section 
VII concludes this paper with a summary of our results. As an application we also compute the 
temperature dependence of the in-medium decay widths to one-loop order for on-shell a and tt 
mesons at rest. 

We use the imaginary-time formalism to compute quantities at finite temperature. Our notation 

is 

lfik)^T f; J^fi2mnT,k) , j J{x) ^ j'J^ dr j £^ f{T,^) . (1) 
We use units fi = c = ks = I- The metric tensor is g^'' — diag(+, — , — , — ). 



II. THE EFFECTIVE POTENTIAL IN THE CORNWALL-JACKIW-TOMBOULIS 

FORMALISM 

The notion of an effective action is quite useful for studying theories with spontaneously broken 
symmetries. For translationally invariant systems, the effective action becomes the effective poten- 
tial. At the classical level, the effective potential is given by the potential energy (density). The 
vacuum (ground) state is given by the minimum of the potential energy. For theories with spon- 
taneously broken symmetry there may exist infinitely many equivalent (degenerate) minima. At 
the quantum level, there are additional terms in the effective potential, corresponding to quantum 
fluctuations. At finite temperature (and finite chemical potential), the minimum of the effective 
potential corresponds to the thermodynamic pressure . 

The common way to compute the effective potential is via the loop expansion . This approach, 
however, becomes problematic for theories with spontaneously broken symmetries. In particular, 
the energy of quasi-particle excitations with small 3-momenta becomes imaginary. The reason is 
that the requirement of convexity for the effective potential is violated. A way to salvage the loop 
expansion is to perform a Maxwell construction which restores the convexity of the effective potential 
pql . Another way to compute the effective potential is via the CJT formalism . As mentioned 
in the introduction, this method resums certain classes of diagrams and has the advantage that the 
energy of the quasi-particle excitations remains real for all values of 3-momentum. 

Consider a scalar field theory with Lagrangian 

£(0) = i9^09'^0-[/(0) , (2) 

for instance, 0"* theory where 

[/(0) = lmV + A04 . (3) 
The generating functional for Green's functions in the presence of sources J, K reads [pO[: 



Z[J,K] = e^I-^'-^l " / ^'^ exp|/[0] -f + 



(4) 



where >V[J, if] is the generating functional for connected Green's functions, = £ is the 
classical action, and 
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) J; 



>K( 



(j){x)J{x) , 
(j){x)K{x,y)(j){y) . 



(5a) 
(5b) 



The expectation values for the one-point function, (/'(x), and the connected two-point function, 
G{x, y), in the presence of sources are given by 



SW[J,K] _ 

'-^4^oix^y)^M] 



(6a) 
(6b) 



One now ehminates J and K in favor of and G via a double Legendre transformation to obtain 
the effective action 



where GK = J^^^ G{x,y) K{y,x). Thus, 

5T[^,G 
S^{x) 
6r\6,G 



J{x) 



K{x,y)(p{y) 



Km \ = -nK{x,y) 
dG[x,y) 2 



(7) 

(8a) 
(8b) 



For vanishing sources, we find the stationarity conditions which determine the expectation value of 
the field (f{x) and the propagator Q{x,y) in the absence of sources: 



5(l){x) 
5T[4>,G] 



5G{x,y) 







. 



(9a) 
(9b) 



Equation (pli| ) corresponds to a Schwinger-Dyson equation for the full (dressed) propagator. It was 
shown in [pof that the effective action r[0, G] is given by 



r[0, G] = m - i Tr (In G'^) - \ Tr {D'^ G - l) + V^\^, G] 
Here, D^^ is the inverse of the tree- level propagator, 



D-'{x,y;cb) = - 



5(j){x) S(j){y) 



(10) 



(11) 



and r2 [0, G] is the sum of all two-particle irreducible diagrams where all lines represent full propa- 
gators G. 

For constant fields 4>{x) = (f), homogeneous systems, and for a Lagrangian of the type given by eq. 
(0), the effective potential V is given hy V — —TT/Q, where fl is the 3-volume of the system, i.e.. 



V[^,G] = U{^) + 1 I lnG-\k) + 



1 



[D-\k;^)G(k)-l]+V2[4>,G] 



(12) 



with 
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D-\k;^) = -k^ + U"{<j>) 
and V2[(j>, G] = —TT2[<j), G]/fl. The stationarity conditions are given by 

SV[4>,G] 



b=ip,G=g 



dV[(t),G] 



. 



SG{k) 

With eq. (|l^), the latter can be written in the form 
where 



SG{k) 



(13) 

(14a) 
(14b) 

(14c) 
(15) 



is the self energy. Equation (14c) is the aforementioned Schwinger-Dyson equation. The thermody- 
namic pressure is then determined by 

p = ~v[ip,g], (16) 

which, in the absence of conserved charges, is (up to a sign) identical to the free energy density. 



III. THE 0{N) MODEL 
Let us now turn to the discussion of the 0{N) model. Its Lagrangian is given by 



(17) 



where (j) is an A^-component scalar field. For H — and > 0, the Lagrangian is invariant under 
0{N) rotations of the fields. For H — and < 0, this symmetry is spontaneously broken down 
to 0{N — 1), with iV — 1 Goldstone bosons (the pions). The phenomenological explicit symmetry 
breaking term, H, is introduced to yield the observed finite masses of the pions. Spontaneous 
symmetry breaking leads to a non-vanishing vacuum expectation value for (p: 



|(0)| =(/)>0 
assumes the role of in section |l^.) At tree level, 



Nm^ 2 

-r-. 7= cos - 

4A V3 3 



9 = arccos 



HN 



12A 



8 A V A^TO^) 



3/2' 



(18) 



(19) 



For H — 0, cos(0/3) = The inverse tree-level sigma and pion propagators are given by 

D-\k;cf>)^-k^+m' + ^^cl,\ (20a) 
D~^(k;cf,) = -k^ + m^ + ^cb^ . (20b) 

This leads to the zero-temperature tree-level masses 
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2 2^12A/,^ 



2 2 I 



N 
TV 



(21a) 
(21b) 



for the sigma meson and the pion. At tree level, the parameters of the Lagrangian are fixed such 
that these masses agree with the observed values of mo- = 600 McV and — 139 MeV. Then, the 
coupling constant is 



A 



iV {ml - ml) 



where = 93 MeV is the pion decay constant and 

n ml -3 ml 



(22) 



(23) 



The explicit symmetry breaking term is H ~ m^/^. These tree-level results may change upon 
renormalization. 

The CJT effective potential for the 0{N) model is obtained from eq. ( p^ as 
F(0, G,, G,) = im202 + „ 



+ i ^ [\nG-\k) + D-\k- 0) G.(fc) - 1] 



N ~ 1 



2 Jk 

V2{4>, Ga, Gtt 



[lnG-\k)+D-\k;^) G,(fc) - l] 



(24) 



where V2{4>, G^, G^) denotes the contribution from two-particle irreducible diagrams. In the follow- 
ing we include only the two-loop diagrams shown in Fig. |^ in V2 . These diagrams have no explicit (j) 
dependence. Then, using eq. (|l^ only tadpole diagrams (with resummed propagators) contribute to 
the self energies. As explained in the introduction, this corresponds to the Hartree approximation. 
The Schwinger-Dyson equations for the full propagators contain no momentum dependence. Thus, 
these equations are simply gap equations for the masses of the sigma meson and pion. 

On the two-loop level there exist, however, two more diagrams, cf. Fig. ^, which will not be taken 
in our analysis. They depend explicitly on (j) and introduce an additional momentum dependence 
in the Schwinger-Dyson equations, which makes their solution more complicated. However, in the 
large- limit these terms are a priori absent, because they are of order 1/A^. 

In the Hartree approximation, 



V2{4', Gcr,G.^) = 3 



N 



1 2 



G4k) 



+ {N + l){N-l) 



GAk) 



+ 2(N - 1) — 



G„{k) 



G^ik) 



(25) 



The coefficients in this equation are chosen such that, when computing the self energies from eq. 
( |l5| ) and replacing the dressed propagators by the tree- level propagators, one obtains the standard 
results for the perturbative one- loop self energies |23[. 
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FIG. 1. The Hartree contributions to the CJT effective potential. Full lines correspond to Ga, while 
dashed lines correspond to Cr- The four-particle vertex ~ A is represented by a full square. 




/ 



\ 



FIG. 2. The neglected two-loop diagrams for the CJT effective potential. The three-particle vertex ~ 
is represented by a full circle. 



IV. THE STATIONARITY CONDITIONS FOR THE EFFECTIVE POTENTIAL 



The stationarity conditions (14a), (p_4b|) read 



4A f 

g-\k) = D-\k;^) + - J [g4q) + {N + i)g^iq)] . 



(26a) 
(26b) 
(26c) 



The integrals on the right-hand side of the last two equations correspond to the sigma meson and 
pion self energies. According to eq. (p^), they originate from the diagrams of Fig. |^ via cutting one 
of the two loops in these diagrams. As one observes, these terms are independent of the momentum 
k^^ appearing the propagator. The only k dependence on the right-hand side enters through 
and D~^, cf. eqs. (20a) and ( pOb| ) . Therefore, one is allowed to make the following ansatz for the 
full propagators: 



5<T,7r(fc) 



1 



(27) 
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where now and are the masses dressed by interaction contributions from the diagrams of 
Fig. 0. Note that the diagrams in Fig. ^ have an exphcit dependence on the external momentum; 
including them would invalidate the ansatz (p7|). 

The dressed sigma and pion masses are then determined by the following gap equations 



Ml = m2 + ^ [3 ^ 3 Q{M^, T) + {N - 1) Q{AU, T)] 



Mi 



4A 



— [^' + Q{A'U,T) + iN 



l)QiAU,T)] 



Here we introduced the function 

g(M,T) ' ^ 



+ M2 



(27r)3 ek(M) iexp[ek(M)/T]-l 2 



(28a) 
(28b) 

(29) 



where ek(M) = \/k^ + AP. The last term in the integral is divergent and requires renormaiization. 
This will be discussed in section The standard practice, however, is to ignore this term, claiming 
it is independent of temperature. This is wrong, because e]^{M) depends on T through the gap 
equation for M. As is shown below, the correct renormaiization procedure changes the results. 
Finally, ip is determined by 



- <P jm' + ^ [<p' + 3 Q{Ma, T) + {N~ 1) Q(M,, T)] | 



Using eq. ( 28a ) this can be written in the compact form 

H = (p 



(30a) 



(30b) 



Note that this equation does not require ex plici t r enor maiiza tion, and therefore is valid independent 
of the renormaiization scheme. Equations ( ^8a| ), (28b), and ( 30b ) are the stationarity conditions in 
the Hartree approximation. In the case where chiral symmetry is not explicitly broken, H = m^r = 0, 
they imply the following: 

1. ip ^ 0. This is the phase where chiral symmetry is spontaneously broken. From eq. ( pOb| ) 
follows 



(31) 



On the other hand, eqs. (28b) and (30a) can be combined to give 

8A 

Ml^ — [Q{AU,T)-Q{AU,T)] . 



(32) 



This implies that Goldstone's theorem cannot be satisfied in the Hartree approximation at all 
temperatures: Mo- 7^ on account of (^l|), therefore = is not a solution of (3^). (Note, 
however, that after proper renormaiization of the function Q(M, T), M^r can be chosen to be 
zero at one particular temperature, for instance T = 0, but then will be nonzero for other 
values of r.) 

2. (/? = 0. In this phase chiral symmetry is restored and eqs. ( ^8a| ) and ( |28b| ) can be combined to 

8A 



Mi - AI. 



/r2 



N 



[Q{AU,T)~Q{AU,T)] , 



(33) 



which has the solution M„ = M^; the masses become degenerate. 
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Let us now turn to the discussion of the large- A'' approximation which is in fact the N ^ 1 hmit 
of the Hartree approximation. To derive the large- TV limit from the previous results, one simply 
neglects all contributions of order 1/N. Note, however, that ip^ ~ N, cf. eq. (^9|). Therefore, in the 
large- iV limit, the stationarity conditions read 

KC = m^ + ^[3^'' + NQ{AU,T)] , (34a) 
M^ = m^ + ^[^^ + NQ{AU,T)] , (34b) 
H = ^ |m2 + ^[^^ + N Q{M^,T)] \ . (34c) 



This leads to 



and 



= A4 + ^ ^^ (34d) 



M^ip = H . (34e) 

These two equations are valid independent of the renormalization scheme. The latter equation 
implies the following in the case that chiral symmetry is not explicitly broken, H = TO^r = 0: 

\. Lp ^ Q. In this phase of spontaneously broken chiral symmetry Af^r has to vanish, i.e.. Gold- 
stone's theorem is respected in the large- limit. A/o- obeys the same relation (|3^) as in the 
Hartree approximation. 

2. ip = 0. In this phase of restored chiral symmetry we again have = M.^, as in the Hartree 
case. 



V. THE RENORMALIZED GAP EQUATIONS 

As mentioned above, the last term in the integrand in (p9| ) is divergent and requires renormaliza- 
tion. In the following we discuss renormalization with a three-dimensional momentum cut-off (CO) 
and via the counter-term (CT) renormalization scheme. 

In the literature one often encounters the argument that this divergent term does not depend on 
temperature and can therefore either be absorbed in the definition of the renormalized vacuum mass 
in the CO scheme, or it is completely cancelled by a counter term in the CT scheme. This argument 
is correct to one-loop order in perturbation theory, since then the mass M in this term is simply the 
bare mass and independent of temperature. However, in a self-consistent approximation scheme, 
like the Hartree approximation, this argument is incorrect, since the mass M is the resummed mass, 
which becomes a function of the temperature through the self-consistent solution of the gap equation. 
Therefore, removing the divergence in either the CO or CT scheme may leave a finite, temperature- 
dependent contribution. Another way of stating this fact is that, as mentioned in the introduction, 
the renormalization constants may have to be chosen such that they depend on properties of the 
medium, like the temperature. 



A. CO scheme 



The simplest way to regularize the divergent integral is to introduce a three-dimensional ultraviolet 
momentum cutoff, A. Computing the divergent integral then proceeds as follows. 
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,(M) 



1 



(27r)3 2ek(M) 4^^ 
In the limit A — > oo, this yields 



dk 



fc2 



1 



ek(M) 87r2 



AeA(M) -M^ In 



2,„ A + eA(M) 



Q(Af,0) = lim Qa{M) 

A^oo 



d^k 1 

(2^)3 2ek(M) 



/i - A/2 /2 + ^ In 



M2 ^ M2 



where we have introduced a renormalization scale, /i, and following ||27| we have defined 

A— >oo iDTT^ /i^ 

The renormalization is carried out by introducing new parameters |^4| 

m2 4(iV + 2) 
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m 


Xr 


^ ~ 


1 


1 







h 



N 



where m% and Xr are the finite, renormalized mass and coupling constant. 



■ (35) 
(36) 

(37a) 
(37b) 

(38a) 
(38b) 



1. Hartree approximation 

In the Hartree approximation, this leads to the following renormalized gap equations for the sigma 
and pion masses: 



Mi 



+ ^ ^ + nM.^ T) + iN- 1) P(M., T)] 



Mi = 



2A 

iVA^ 

R + 

2A 



TO 



AT A 



R 



N N 
Ml - ml 
4Ak N + 2 

Ml - ml 



4Afl 
TV 



(A^ + 2) [(^2 + p(M,^r)] 



[(^2 + P(M,, T) + {N- 1) P(M,, T)] 
4Afl 



N 



{N + 2)P{M,,T) 



where the function P{M, T) is defined as 
P{M, T) 



M2 , M2 
In — TT 



167r2 



1 



1 



(2^)3 ek(M) exp[ek(Af)/r] - 1 



(39a) 



(39b) 



(40) 



Equations (39a) and ( |39b| ) are equivalent to eqs. (13), (14) in |lj] after the replacements A 
A/6, 02 ^ ^ p{M,T) P/[M]. [Note that the terms - Ml ~ Ml in eqs. (13), (14) of (l| can 
be eliminated by taking the difference of eqs. (13) and (14).] 

In the limit A — > cxd, A — > 0^ in order to have a finite A^j, and the (bare) theory becomes unstable 
(see also Ref. [^). Also, the (renormalized) masses obey Ml — Ml, cf. ( ^9aD , ( |39h| ), which is 
undesirable. It would imply that chiral symmetry is unbroken, even when (p ^ 0. This problem was 
also addressed by the authors of ||T^. On the other hand, for < A < cxd, A^j ^ O"*" in the limit 
A — > 0, indicating that the (renormalized) theory becomes trivial ]25t . 

Therefore, in the CO scheme the gap equ ati ons c an only be meaningfully studied for finite A. In 
this case, in the original gap equations (28a), (28b) we replace 
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Q(M, T) ^ Qa(M) + Qt{M) , (41) 

where 

/(Ph 1 1 
{2TTy ek(M) exp[ek(M)/TJ - 1 

The integral Qt{M) is UV- finite and does not require the introduction of a momentum cut-off. 
Consequently, the gap equations read 

Ml = ^ ^ {3(^2 ^ 3 [Qt{M„) + QA(Af,)] + (iV - 1) [Qt{M^) + QA(Af,)]} , (43a) 

Ml =m^ + ^{ip^ + [QAM,) + Qa(M,)] + (iV + 1) [Qt(M,) + Qa(M.)]} . (43b) 

The cut-off A has to be determined from the values of and Af^ at T = 0: 

4A 

=m2 + — [3/2 + 3QA(m,)-f (A^-l)gA(m,)] , (44a) 
^1 = 171^ + ^[fl + Q!y{m„) + {N + l)QA{m^)\ , (44b) 



wher e we have used Lp{T ~ 0) — fj,- In the chiral limit {H ~ m^r — 0), the difference of (44a) and 
(HI) reads 

8A 

= ]^ [/' + Qa(^"-) - Qa(0)] • (45) 

However, from the stationarity condition ( t30b| ) we conclude that = 8\fl/N in the chiral limit. 
This immediately leads to 

QA(m,) = gA(0) , (46) 

which for finite m„ can only be fulfilled if A = 0. This, however, is exactly the case treated in pl| , 
without renormalization. 

In conclusion, the CO scheme fails to provide a consistent renormalization of infinities in the phase 
of broken chiral symmetry in the Hartree approximation when H — to^ = 0. Note that the same 
conclusion can be reached with a four-dimensional momentum cut-off. This failure can be traced to 
the fact that in the Hartree approximation diagrams of the type shown in Fig. || are not included 
(cf . the perturbative renormalization of the linear sigma model , see also the discussion in ) . 

Due to this failure, no results will be shown for the Hartree approximation with CO renormaliza- 
tion. However, we note that the case of explicitly broken symmetry, H ^ 0, tUt^ > 0, is free of this 
problem. Then, the difference of eqs. (44a) and (44b) determines the coupling constant as 

A^A(A)^^ ^. — ^ ^ . (47) 

8 + QA[ma) - (3A(TO7r) 

The mass parameter is given by 

= - ^ (^ + 2) QaM ■ (48) 

H is determined from (^) to be iJ = f^[ml ~ 8X{A) fl/N]. 
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2. Large-N approximation 



In the large- iV limit, the renormalized gap equation for the pion mass reads 



Ml = ml 



(49) 



while is still given by ( 34d ). This again has the consequence that Ala- — in the limit A ^ oo, 
i.e., A — > 0^, which as discussed above is an unwanted feature. On the other hand, there is no 
inconsistency in t he la rge- iV approximation for finite A. A is a free parameter and the gap equations 
to be solved are ( 34d ) for the sigma mass and 



M2 ^ m'' + ^ {if'' + N [Qt{AU) + Qa{AU)]} 



(50) 



for the pion mass. The parameters are again determined from M„{T = 0) = nia, Mt^{T = 0) = m^r, 
and (/)(T — 0) — /tt- From these conditions we derive that the coupling constant is still given by its 
tree-level value, eq. (E2|), but is now determined from 



3 m 



^-4AQa(to,) . 



(51) 



H retains its tree-level value on account of (34e). 



B. CT scheme 



In the CT scheme, counter terms are introduced to subtract the UV divergences in Q(M,T). 
Rewrite 



1 



(27r)3 2ek(M) J {2tt)^ + 



1 



(52) 



where /co G R- with k'^ ~ k^ + k^, d^/c = d^lidko. To determine the counter terms, expand the 
integrand in a Taylor series around ~ fp, where /i is the renormalization scale. 



1 



fc2 + M2 fc2 + ^: 



- M= 



fc2 



(53) 



The d'^k integral over the n — term in this expansion is quadratically divergent, while the integral 
over the n — I term diverges logarithmically. The counter terms are chosen to remove these two 
terms, such that the renormalized result for the divergent integral is 



d'^k 
(2^)4 



1 



1 



M 



2 1 



Af2 



+ ^2 (-^2 ^^2)2 



oo „ 

= 5:(M^-M^)"/ 

n=2 



1 



(27r)4 (fc2-hAi2) 



2\n+l 



(54) 



Note that the second counter term depends on the temperature through M. This fact represents 
the aforementioned possibility of having temperature-dependent counter terms in self-consistent 
approximation schemes, and was already discussed by the authors of ||l0|] . They also pointed out 
that this problem does not occur in less than three spatial dimensions. This is obvious from eq. 
(p^), because then the second counter term is finite, and thus not required. In contrast, cither in 
ordinary perturbation theory or in optimized perturbation theory [ p7[ renormalization at T = is 
sufficient to remove all divergences. 

The last integral in (p3) is finite and equal to 
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1 



1 ^, 



2(l-n) 



Expression (EJ) can be rearranged to give the final result 



_<£k_ 
(2^ 



//-2 



1 1 ^i^ -M 



(4^) 



Af2 In ^ - + ^2 



To obtain the renormalized gap equations, simply replace Q{M,T) as given in (29) by 

QiM,T)=QTiM)+Qf,{M) , 



where 



{Any 



A/2 In 



M2 

7^ 



(55) 



(56) 



(57) 



(58) 



The renormalization scale is chosen to give the correct values for sigma and pion mass at T = 0. 

As an alternative to the above procedure, one can also compute ( p2| ) in dimensional regularization, 
i.e., in d space-time dimensions, where the coupling constant g is replaced by gfl'^. Here, fl is the 
renormalization scale in dimensional regularization and e = 4 — d. In order to obtain (p6|), one has to 
add a counter term M'^/{8TT'^e) + ij.'^ / {16tt'^) . Here, /i is the renormalization scale from the previous 
treatment and related to /i by /i^ = ATTe'^js^, where 7 is the Euler-Mascheroni constant. Note again, 
that the counter term depends implicitly on the temperature through the resummed mass M . In 
the Appendix, we furthermore show that the CO and CT schemes are equivalent for unbroken 0{N) 
symmetry. 



1. Hartree approximation 
In the Hartree approximation, the gap equations read 

Ml ^ + ^ {3 ^2 _,_ 3 [Q^(M^) + Q^(A/,)] + (iV - 1) [Qt{M^) + Qt,{M^)]) , (59a) 
A/2 = m^ + ^{^^ + [QriM,) + Qt.{M„)] + (iV + 1) [Qt(A/.) + Q^(A/,)]} . (59b) 
The renormalization scale ^ is determined from the vacuum values for the sigma and pion masses: 



ml^m^ + ^ [3/2 + 3 Q^(to.) + (N - 1) Q^(m,)] 
ml = m^ + ^ [/2 + Q^(m,) + (A^ + 1) Q^(to.)] . 



In the chiral limit, the difference of these two equations reads 

f 



8A 

TV 



167r2 fi^e 



(60a) 
(60b) 

(61) 



However, in order to be consistent with the (generally valid) equation ( pOb| ) , there is only a single 
choice for the renormalization scale, jj.^ = m1/e. Then, the coupling constant is given by its classical 
value, A = Nml/{S,f^), while 



2 N ^ ' 167r2 



(62) 
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In the case that chiral symmetry is exphcitly broken, the difference of ( |60a| ) and (60b) yields the 
following equation for the coupling constant: 



8 + [m^ hi{m'^/ ijpe) — \n{m'^ / fj? e)] /IBtt^ 



A(/i) 



(63) 



i.e., A runs with the renormalization scale. However, there is one value for the renormalization scale, 
where A retains its tree-level (i.e. classical) value, 



= Mci = cxp 



(Inrn^ — 1 ) 



(Inm^ - l) 



(64) 



The results for the Hartree case with explicitly broken symmetry presented in the next section will 
exclusively employ this value of /x. The mass parameter is determined from 



ml- 3 ml 



4A 



(65) 



H can be obtained from (34e) at T = 0. 



2. Large-N approximation 
In the large- iV limit, the gap equations to be solved are ( |34d| ) for the sigma meson and 

= + ^ { (^2 ^ TV [Qt {M^ ) + {M^)] ] (66) 

for the pion. In this case, /i is a free parameter, and cannot be fixed by the vacuum values for the 
sigma and pion masses. A and H are always given by their tree-level values. The mass parameter is 
determined from 

m2 = -^^i-^-4Ag^(mO . (67) 



VI. RESULTS 



In this section, we discuss numerical solutions of the gap equations for the meson masses and the 
stationarity condition on Lp. Three different cases are considered: the large- approximation in (a) 
the CO scheme, (b) the CT scheme, and (c) the Hartree approximation in the CT scheme. The 
Hartree approximation in the CO scheme will not be discussed, due to the problems exhibited in 
section We focus separately on the cases TOtt = and > 0. 



A. = 



Figures || (a,c,e) show the meson masses and (b,d,f) </? as functions of temperature. Results 
for the large- iV approximation with CO renormalization are shown in parts (a,b), and with CT 
renormalization in (c,d). Results for the Hartree approximation with CT renormalization are shown 
in (e,f). For comparison, the dashed lines in each figure correspond to the unrenormalized results of 

In Figs, (a,b), in the phase of spontaneously broken symmetry, there is no difference between 
the unrenormalized and renormalized cases. To understand this, first remember that = 0, cf. the 
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discussion following eq. (34e ). T herefore, on account of (B4d), Mg. is simply given by {SX/NY^'^tp. 
In turn, (p is determined by (34c). However, for Mj^ — 0, this has the simple form 







4A 

iV 



4A 



12 



Qa(0) 



Using ( ^l|) for jti^ — 0, this becomes 







to; 
2 



4A 



4A; 



y2 



(68) 



(69) 



N ^ 12 ' 

which is the same condition as in the unrenornialized case (where Qa is absent). Since the coupling 
constant is given by its tree- level value (|2^), this immediately leads to the conclusion that the 
temperature for chiral symmetry restoration is 



(70) 



In the restored phase, = 0, sigma and pion masses are equal, and given by eq. (34a) or (34b). 
These equations are cut-off dependent, on account of (^). The mass is decreasing for increasing A. 

In Figs. H (c,d), we show results for the large- TV approximation in the CT scheme. In the broken 
phase. If > 0, renormalization again does not affect the masses or ip. In the phase of restored 
symmetry, p = 0, the sigma and pion masses are degenerate, but depend on the renormalization 
scale. They decrease for increasing ^. Note the similarity between the masses in the CO and the CT 
scheme when choosing the same value for the cut-off A and the renormalization scale /x. Considering 
that both renormalization schemes are fundamentally different, this similarity is quite surprising. 
Another important conclusion is that renormalization of the gap equations does not destroy the 
second-order nature of the transition. 



Large N, CO Large N, CT Hartree, CT 




100 200 300 100 200 300 100 200 300 400 



T [MeV] T [MeV] T [MeV] 

FIG. 3. The meson masses and </5 as a function of T for m-n- = 0. 
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The results in the Hartree approximation, eqs. (28a) - ( ^Ob| ), are displayed in Figs. ^ (e,f)- As 
in the unrenormalized case, we obtain a first order transition, with a transition temperature that 
appears to be slightly higher than in the unrenormalized case. To determine this temperature, 
however, one would have to analyze the shape of the effective potential, which is outside the scope 
of this paper. 



B. = 139 MeV 

Fig. ^ (a,c,e) shows the temperature dependence of the meson masses and Fig. I (b,d,f) the 
function (p{T) in the case of explicit symmetry breaking. As in Fig. |3[ large- iV results are shown in 
(a,b) for the CO scheme and in (c,d) for the CT scheme. Parts (e,f) show our results for the Hartree 
approximation with CT renormalization. As already observed in the chiral limit, there is a striking 
similarity between the results in the CO and the CT schemes when choosing A = fi. Also, increasing 
the cut-off or the renormalization scale tends to increase the temperature at which (approximate) 
symmetry restoration takes place. 

Baym and Grinstein [p^ noted that the additional terms originating from renormalization have 
the effect that the gap equations do not have a solution beyond a certain temperature (see also 
1^,^). We found evidence for this in the CT scheme at temperatures above 400 MeV. In the CO 
scheme with a finite A, this phenomenon does not occur. 
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FIG. 4. The meson masses and <^ as a function of T for =139 MeV. 
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VII. CONCLUSIONS AND OUTLOOK 



In this paper, we have studied the temperature dependence of sigma and pion masses in the 
framework of the 0{N) model. The CornwaU-Jackiw-Tombouhs formahsm was apphed to derive 
gap equations for the masses in the Hartree and large-iV approximations. Renormahzation of the 
gap equations was carried out within the cut-off and counter-term renormahzation schemes. In 
agreement with jl^, it was found that the cut-off scheme is flawed when the cut-off A ^ oo. We 
therefore studied this renormahzation scheme for A < oo. For the Hartree approximation we found 
that, in the chiral limit (171^^ ~ 0), there is no finite value for the cut-off, which is consistent with 
the set of stationarity conditions for the effective potential; A = is the only possible choice. This 
problem was not encountered in the large- iV approximation; here any choice for A is possible. In the 
counter-term renormahzation scheme, the Hartree approximation can be consistently renormalized, 
but in the chiral limit, the renormahzation scale is restricted to a unique value in order to achieve 
consistency with the stationarity conditions for the effective potential. In the large-iV limit, the 
renormahzation scale can be chosen arbitrarily. Changing the cut-off in the cut-off scheme or the 
renormahzation scale in the counter-term scheme changes the meson masses at a given temperature. 
The reason is that, in the self-consistent approximation schemes considered here, the renormahzation 
constants (or counter terms, respectively) may depend implicitly on temperature. This does not 
occur when renormalizing ordinary perturbation theory. 

Our results can be compared to those of Roh and Matsui ||l^ and Chiku and Hatsuda [0 . The 
authors of ]l5| computed the sigma and pion masses from the second derivative of an effective 
potential which was determined via the standard loop expansion approach. Being aware that this 
approach fails for theories with spontaneously broken symmetry, they corrected the resulting expres- 
sions to obtain gap equations which look similar to the ones in the Hartree approximation. (They are 
identical to Baym and Grinstein's modified Hartree approximation pO| ] ) . The stationarity condition 
for if, however, was taken to be the same as in the large-A^ approximation. Thus, their solutions 
respect Goldstone's theorem in the phase of spontaneously broken symmetry, similar to the large- A'^ 
approximation discussed here, while the transition in their model is first order (in the chiral limit), 
like in the Hartree approximation. 

The authors of employ optimized perturbation theory to compute the sigma and pion masses. 
This approach has the advantage that renormahzation is straightforward. The results are similar to 
those of 

Recent dilepton experiments at CERN-SPS energies |^ have generated interest in medium mod- 
ifications of meson properties such as their mass and decay width. In general, the meson mass 
(squared) is given by the inverse propagator at = 0, = G^^{Q)- The decay width of a particle 
with energy lu at rest, k = 0, is given by ^{uj) = — ImH(a-', 0)/a; |29[ , where n(cj,k) is the self 
energy. In the CJT formalism, Q^^{k) = D^^{k; ip) + H(fc), cf. eq. ( |14cD . In the Hartree or large-A^ 
approximation studied here, the self energies do not acquire an imaginary part, because they are 
simply constants, and thus only shift the mass of the particles. Therefore, in these approximations, 
the particles are true quasi-particles with vanishing decay width. This would change if we included 
the diagrams of Fig. || in the effective potential, because, as is well-known |29|, the imaginary part 
of these diagrams corresponds to decay and scattering processes. 

To include these diagrams in the above treatment, however, is prohibitively difficult, because 
then the simple momentum dependence of the propagators Gcr,-K{k) in eq. (27) changes, since the 
self energies become explicitly momentum dep ende nt. Then, instead of simple gap equations for 
the meson masses, the stationarity condition ( |l4b| ) becomes an (infinite) set of coupled integral 
equations for the propagators Qa,-K{k). 

Therefore, as a first approximation, we compute the decay widths from the self energies corre- 
sponding to these diagrams, but with internal lines given by the Hartree or large- A'^ propagators 
( p7| ) . This is equivalent to computing the decay width to one-loop order in perturbation theory, but 
taking the medium-modified masses of the particles computed above instead of the vacuum masses. 
The on-shell decay width of a and tt mesons at rest is then given by the following expressions , 
valid for 2M^<M„: 
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7<T = — ^ 1 1 f coth — , (71a) 



_ /4A0y / 4M2 l-exp[-MjT] 1 

^ TV y STTMsy M2 1 - exphAf2/2m,r] exp[(A/2 - 2Af2)/2M,r] - 1 ■ ^ > 

These quantities are shown in Fig. ^ for ttItt — 0, and in Fig. ^ for = 139 MeV, for the cases 
discussed in Figs. ^ and ^. For = and in the large- iV approximation, pions are true Goldstone 
bosons, and therefore their decay width vanishes befow the temperature corresponding to chiral 
symmetry restoration, see Figs. |^ (b,d). This is different in the Hartree approximation, where 
Goldstone's theorem is violated, cf. Fig. ^ (f), and when chiral symmetry is explicitly broken, Figs. 
^ (b,d,f). The reason is that, because pions have a finite mass, they can acquire a finite decay width 
on account of the absorption processes ttct tt and tttt — » a. For massless particles, these processes 
are kinematically forbidden. 

Sigma mesons, however, can always decay into two pions, and therefore acquire a large decay 
width, cf. Figs. ^ and ^ (a,c,e). All decay widths vanish above the temperature where Ma- becomes 
smaller than 2 A/jr- This, however, is an artefact of the one-loop approximation. In two-loop order, 
the scattering processes aa — > aa, aa — s- tttt, cttt — s- crvr, and tttt tttt lead to a finite decay width 
for all particles even above this threshold. 

The decay widths and masses computed here are relevant for the formation of disoriented chiral 
condensates [Q, since they enter the evolution equations of the long- wavelength modes. This will 
be the subject of a subsequent investigation js^. 
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APPENDIX A: SCHEME EQUIVALENCE FOR UNBROKEN 0{N) SYMMETRY 



In this Appendix, we show that the CO and the CT schemes are equivalent when the 0{N) 
symmetry is not broken (ip = 0). In this case, the gap equations become degenerate. 



M2 = m^ + ^{N + 2) Q{M, T) . 



(Al) 



(We only discuss the Hartree case here, for the large- iV approximation, simply replace TV + 2 by iV.) 
In the CO scheme, using eqs. p9) and (pq), this becomes 



4Aj 
N 



{N + 2) 



QAM) 



In — 7T- 



(A2) 



The renormalization scale /i^ can be determined from the T — limit of this equation. For unbroken 
0{N) symmetry, M(T = 0) = rriR, which then yields — ttir. 
On the other hand, in the CT scheme, we have 



2 , 4Afl, 

mR + —{N- 



Af2 ^ AP 
In — 7T- 



(A3) 



Here, we made the finiteness of m and A explicit by replacing them with mji and A^. Again, the 
condition M{T = ) = mn yields /i ~ mn. 

The last term in ( |A3[ ) leads to an apparent difference between the two schemes. However, shifting 
the coupling constant by a finite amount. 



1 1 



4{N + 2) 



(A4) 



one obtains (A2), which proves the equivalence of both schemes a fter properly redefining the coupling 
constant. The same conclusion can be reached starting from (Al) and using instead of ( |3^ ) the 
modified renormalization conditions 



1 

A^ 



A{N - 



1 4:{N + 2 

^ 167r2iV 



m2 4(iV + 2) 
X ^ N 
1 , 4(iV + 2) ^ 
A 



h , 



N 



(A5a) 
(A5b) 



which then leads to (A3). 
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